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1 Introduction 



The understanding of the nonperturbative aspects responsible for the long distance be- 
havior of confining Yang- Mills theories is one of the most challenging issues in theoretical 
physics. 

The common paradigm in the quantization of non-Abelian gauge theories relies on the 
Faddeev-Popov formula, which enables us to construct a meaningful version of the theory 
at the quantum level in the weak coupling regime, i.e. a theory which is renormalizable and 
unitary. We underline that unitarity is here referred to the case in which asymptotic fields 
can be safely introduced, so that a physical asymptotic subspace with positive norm can 
be constructed through the use of the BRST symmetry. The situation becomes very differ- 
ent when dealing with the infrared behavior of a confining theory as, for example, QCD, 
where asymptotic quark and gluon fields cannot be introduced due to color confinement. 
Despite the progress done in the last decades, a complete understanding of the behavior of 
non-Abelian gauge theories in the infrared regime accounting for color confinement is still 
lacking. 

One possible way to face this problem is to go back to the Faddeev-Popov quantization 
scheme and realize that, as it stands, it does not account for the existence of equivalent 
gauge field configurations fulfilling the same gauge fixing condition. This is the Gribov 
problem pQ. Let us quote, as an example, the case of the Landau gauge, which has been 
intensively investigated in the last years. We remind that, according to Singer's result [2], 
the existence of the Gribov phenomenon is a general feature of the gauge fixing procedure. 
Gribov copies have in fact been observed in most of the employed covariant and renormal- 
izable gauges, namely: the Landau gauge [H EJ EJ EJ [7] , the class of general linear covariant 
gauges [8], the maximal Abelian gauges [U HQ]. Also, the Coulomb gauge provides an 
example of a non-covariant gauge exhibiting the Gribov phenomenon \12\ 113] . The 
fact that the Faddeev-Popov formula is plagued by the existence of Gribov copies means 
that the gauge fixing procedure does not pick up a unique representative for each gauge 
orbit. As a consequence, the counting of the relevant field configurations is not properly 
done. Although a complete resolution of the Gribov issue is not yet available, it has be- 
come now clear that taking into account the Gribov copies requires a modification of the 
Faddeev-Popov quantization formula in the infrared non-perturbative region, as expressed 
by the so called Gribov- Zwanziger action |14} I15j. which enables us to restrict the domain 
of integration in the functional integral to the region in field space bounded by the first 
Gribov horizon. This restriction deeply modifies the behavior of the Green's functions 
of the theory in the infrared. As an example, the two-point gluon correlation function 
turns out to be suppressed in the infrared, exhibiting complex poles. As such, it cannot 
describe physical excitations, i.e. gluons are destabilized by the presence of the Gribov 
horizon and are not part of the physical spectrum, as expected in a confining Yang-Mills 
theory. Recently, a refined version of the Gribov- Zwanziger theory has been established 
[EJ CE7J HB], through the introduction of dimension two condensates. We point out that 
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the predictions of the refined version of the Gribov-Zwanziger theory turn out to be in 
remarkable agreement with lattice numerical simulations, both at the level of the gluon 
and ghost propagators, see |19l [20] and references therein, as well as at the level of the 
spectrum of the lightest glueball states [21] . All these nontrivial results can be taken as 
evidence that the Gribov phenomenon plays a rather relevant role in the infrared physics 
of Yang-Mills theories. 

The aim of this work is that of pursuing the study of the issue of the Gribov copies, 
already started in [10]. More specifically, we shall provide examples of infinite classes of 
renormalizable zero modes of the Faddeev- Popov operator in the Landau gauge, d{Af = 0, 
i = 1,2, ..,d, in d = 2,3,4 Euclidean dimensions. We shall look thus at normalizable 
solutions of the equation 

M ab u b = , (1) 
where A4 ab stands for the Faddeev-Popov operator in d dimensions 

M ab = -(d 2 8 ab +e acb A c i d i ^J , (2) 

and e acb are the structure constants of SU(2). Nontrivial solutions of the equation (fTJ) im- 
mediately lead to the existence of the Gribov copies. Considering in fact the two equivalent 
gauge configurations (Af,Af) 

Af = Af- (6 ab di + e acb A^j u b , (3) 

it turns out that, due to eq. ([!]), both fields Af and A? fulfill the Landau gauge condition, 
i.e. 

dill = diAf + M ab uj b = diAf = . (4) 
Therefore, the gauge configuration Af is a Gribov copy of Af. 

The strategy which we shall follow to obtain nontrivial normalizable solutions of eq.([T]) 
is that outlined by Henyey in [3], and already successfully employed in [ID]. Henyey's 
strategy amounts to find a parametrization of both zero mode uj a and related field configu- 
ration Af in such a way that eq.([I]) becomes an algebraic equation enabling one to express 
the gauge field in terms of the function parametrizing the zero mode itself. Despite its 
simplicity, Henyey's framework turns out to be very powerful in order to provide examples 
of gauge copies with the desired asymptotic behavior for the resulting gauge field configu- 
ration. 

In what follows we shall consider two kinds of boundary conditions for the behavior of 
the gauge field at infinity. The first type is expressed by the requirement that the Yang- 
Mills action is finite, i.e. 

Sym = ^2 j ddx F ii F ?i < 00 • ( 5 ) 
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The finite action condition looks a natural requirement to be imposed, given the impor- 
tant fact that the action has the meaning of Boltzmann weight in the Euclidean functional 
integral. 



The second type of boundary condition which we shall discuss amounts to require that 
the behavior of the gauge field at infinity is such that the Hilbert norm \\A\\ 2 is finite @], 
namely 



Condition ([6]) is stronger than the requirement of finite action, eq.©. Field configurations 
fulfilling the condition of finite norm | \A\ | 2 automatically lead to a finite action, the converse 
being not true, as it will be shown by explicit examples. The relevance of the Hilbert norm 
||^4|| 2 relies on the fact that, in the Landau gauge, several properties of both Gribov region 
and fundamental modular region can be established by demanding that the functional 



is well defined and differentiable along the gauge orbit of A{ |22[ [23"1 [24"] . It turns out in 
fact that one possible way to introduce the Gribov region is to identify it with the set of 
all relative minima of the functional /a [it] in the space of the gauge orbits, while the fun- 
damental modular region corresponds to the set of all absolute minima. The fundamental 
modular region is contained within the Gribov region and enjoys the property of being free 
from Gribov copies [221 1231 124] . Though, a practical way of implementing the restriction 
to the fundamental modular region in the domain of integration in the functional integral 
is, at present, not known. 

The present work is organized as follows. Sections 2,3,4 are devoted to construct explicit 
examples of infinite classes of normalizable zero modes of the Faddeev-Popov operator, 
whose corresponding gauge field configurations display an asymptotic behavior compati- 
ble with the requirements ([5]), ([6]) in d = 2,3,4 dimensions, respectively. In Sect. 5 a few 
concluding remarks are collected. 

J We remind that the Gribov region is defined as the set of field configurations fulfilling the Landau 
gauge condition and for which the Faddeev-Popov operator M. is strictly positive, namely 




(6) 




u'AiU + v)diU 



(7) 



n = { At ■ d l A a l = , M 



ab 



(d 2 5 ab + e acb A1d^ > } 



(8) 
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2 Henyey's construction of zero modes of the Faddeev-Popov 
operator in the Landau Gauge in d = 2 dimensions 

Let us begin by discussing the construction of an infinite class of normalizable zero modes of 
the Faddeev-Popov operator in d = 2 dimensions. Let Aj, i = 1,2, be a SU(2) Lie-algebra 
valued gauge field in Euclidean space, 



A = A^-j , [a a ,a°]=2ie Mc a c , (9) 
where cr a , a = 1, 2, 3, denote the Pauli matrices. 

Constructing normalizable zero modes of the Faddeev-Popov operator amounts to find 
solutions of the equation 

d 2 oo a + e acb A'rd i uj b = , (10) 

with 

J d 2 xu a uj a < oo . (11) 

Adopting polar coordinates (r,</?jl, the Landau gauge condition, diA? = 0, can be fulfilled 
by setting 

A] = 6 b3 ai = 5 b3 e ljXj f(r) , (12) 

so that eq. ffTUj) becomes 

d 2 u a + /(r)e 3ac e ij x i a i w c = . (13) 

The quantity f(r) in eq. (!12p stands for an arbitrary function depending only on the radial 
coordinate r. Looking at the second term, one recognizes the angular momentum L 

L = eijxidj = g^- (14) 

Then, choosing the components of oj a (r, if) along the directions orthogonal to Af in color 
space, namely u a = (a/ 1 ), a/ 2 ), 0), we obtain 

d 2 JV+f(r)LujW = 

02 w (2) _ /( r )LcjW = o. (15) 
We parametrize (w^^u/ 2 )) in the following way 

w^^r, ip) = a(r) sin tp , 

u( 2 \r, <p) = a{r) cos (p. (16) 



2 see (Appendix [SJ 
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where cr(r) is left free. As a consequence, eqs. (fl~6j) merge into a single equation, which 
reads 

a"(r) + -a'(r) - \a(r) = f(r)a(r) . (17) 

Expression (|17p illustrates in a simple way Henyey's construction [3]. Both gauge field, 
eq,(|12p. and zero mode, eq. (|16j) . are parametrized in such a way that eq. (|17f) allows us to 
express the function f(r) in terms of the quantity cr(r). More precisely, making the change 
of variables 

<r(r) = rip{r) , (18) 

it follows that 



■0(r) \ r J 

What remains to be done is to find a suitable function ip{r) which provides a normalizable 
zero mode u a and a corresponding gauge field A°" with the desired asymptotic behavior. 
To that purpose we shall employ the following trial expression for ijj{r): 

Mr) = , (20) 

where the exponent m is positive and a is an arbitrary positive constant. From expression 
(|19p . for the function f(r) one obtains 

f(r) = - 1 ——(r 2m ' 2 (m 2 -2m)-ar m - 2 (m 2 + 2m)) . (21) 

(r m + a.) 1 ' 

Also, the norm of the zero mode u> a is easily evaluated and turns out to be given by 

.(■00 Jl 

72^ , ,a, a 



llwlr = diwVa / cirr- ^ . (22) 

" " J Jo (r™ + a) 2 y ' 

This expression will be convergent for m > 2. Note that in this case we have, asymptoti- 
cally, /(r) ~ l/?" 2 , and thus 

Ai ~ 1/r , /or r^co. (23) 



It is easy to check that the asymptotic behavior (|23p leads to a finite Yang-Mills action. In 
fact, one observes that, since only the 3-color component of the field is nonvanishing, 
eq. (fT2j) . the field strength, Ff- = diA a -—djA a i +e ahc A\A c -, has a unique nontrivial component, 
namely 

4 = d % A) - djAl . (24) 
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Therefore, making use of eq. (|12p . we have 



4 = 2e ij f(r) + e lk x k xi^± - e ]k x k x if -^ . (25) 
Moreover, using the identity 

^kx j + e kj x l + ejiX k = , (26) 

we find 

4 = e y (2/(r) + r/'(r)). (27) 
Finally, for the Yang-Mills action one gets 

Sym = / ^xFf^ = ^1^ drr(2f(r) + r/'(r)) 2 < oo , (28) 

which is obviously finite, as f(r) ~ 1/r 2 . Observe, however, that this is not the case of 
the Hilbert norm ||^4|| 2 , which is logarithmic divergent, due to Ai ~ 1/r. We see thus that 
expression (|20j) gives rise to an infinite class of normalizable zero modes of the Faddeev- 
Popov operator, parametrized by the exponent m > 2, whose corresponding gauge field 
configurations yield a finite Yang-Mills action. 



3 Zero modes in d = 3 

In order to construct examples of zero modes in d = 3, we shall proceed by generalizing the 
set up outlined in the case d = 2. Adopting polar coordinates (r,6,ip), as starting gauge 
configuration A^,i = 1,2,3, fulfilling the Landau condition, diA? = 0, we shall employ 

AI = e ctjXj h(r) , (29) 

where the function h(r) depends only on the radial coordinate r. The differential equation 
for the zero modes reads 

d 2 u a - h(r)e acb e cij x i d j uj b = . (30) 

As already observed in the previous section, one easily recognizes in the second term of 
eq. (|30p the angular momentum L a 

L a = e aij Xidj . (31) 

Moreover, similarly to the d = 2 case, we choose uj a (r,ip,9) = (u/ 1 ), oj( 2 \ 0), so that eq. (|3U|) 
gives 

dV 1 ) - h(r)L 3 ojW = o, 

dV 2) + h(r)L 3 u (1) = 0, (32) 
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where 

L ' = ^- < 33 > 

Therefore, setting 

u)^\r, tp, 9) = <r(r) sin ip sin 9 , 

cj( 2 \r,(p,9) = a{r) cos </? sin 9 , (34) 
it follows that eqs. (|3~2j) reduce to a single equation, given by 

a"{r) + -</(r) - -|<r(r) = h(r)a(r) . (35) 

Adopting thus the change of variables a(r) = rip(r), we obtain 

h(r) = — L (V(r) + V(0 ) . (36) 



?/>( r ) \ r 

In order to find an infinite class of normalizable zero mode we proceed as before and we 
make use of the function 

^ = ^r— > ( 37 ) 

rP + a 



with p an arbitrary positive exponent and a a positive constant. From eq. (|36p we get 

h ^ = {r pl a) 2 (r 2p_ V " 3p) - or*- V + 3p)) . (38) 

The norm of the zero mode w a is easily computed and yields 

roo r 4 

IM| 2 oc / dr- -77 . (39) 

" Jo (rP + a) 2 v ' 

The integral of equation (|39|) is convergent for p > 5/2. Therefore, expressions (|34p and 
(|37p provide an infinite class of normalizable zero modes, parametrized by the values of the 
exponent p > 5/2. 

Also, for a generic value of p > 5/2, the asymptotic behavior of h(r) in expression (|38p is 
/i(r) ~ 1/r 2 for r — > 00. As a consequence, the behavior of the corresponding gauge field 
A1, eq.® , is 

Ai ~ - , for r -> 00 . (40) 
r 

Again, similarly to the previous case in d = 2, this gauge configuration (|38[) gives rise to a 
finite Yang-Mills action. In fact for the field strength we get 

Ffj = 2e aji h(r) - e aik x k x j ^-^- + e ajk x k x i ^-^- + e abc e bik e cj ix k x l h(r) 2 . (41) 
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Making use of the algebraic identity 

expression f)41[) becomes 

F% = e aji (2h(r) + rh'(r)) + ey fc (h{r) 2 + ^\ x a x k . (43) 
Therefore, for the Yang-Mills action one finds 



/ S XF ^ = 4. / ^ [6 (2,(0 + r h ' { r)f - * (2,(0 + rtf (r)) + ^ 



+2r 4 f^ + ,(r) 2N2 



(44) 



which is easily seen to be finite for p > 5/2. 



However, for a generic value of p > 5/2, the asymptotic behavior (I40p does not guarantee 
that the Hilbert norm \\A\\ 2 = f d^xAfAf is finite too. Nevertheless, for the particular 
value of p = 3, the expression for h(r) reduces to 

i 18qt 

so that 

^| p =3"' ^ or r^°o, (46) 



yielding a finite Hilbert norm 
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lp=3 



| d 3 x(^)| p=3 <oo. (47) 



Evidently, the field configuration (j45|) produces a finite Yang-Mills action. Although the 
expression 

= -IT— > ( 48 ) 

gives only a particular example of a normalizable zero mode whose corresponding gauge 
field, eq.(|46|). displays a finite norm ||A|| 2 , it is relatively easy to construct now an infi- 
nite class of zero modes giving rise to a finite norm. To that purpose, we introduce the 
generalized function 

^ = 3^Rn^ > ( 49 ) 

r 6 + pr n + a 
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where f3 is a positive constant and n a positive exponent such that n < 3. From eq. ([36|) . 
for the function h(r) we obtain now 

h(r) = (ra+ ^n +a)a [/3 2 r 2 «- 2 (n 2 - 3n) + /^^(Qn - n 2 - 18) 

-(3ar n - 2 (n 2 + 3n) - 18ra] . (50) 

We see that we must constrain n > 2 in order to avoid a singularity at the origin. Asymp- 
totically we have h(r) ~ \j r ^>- n anc i thus vlj ~ l/r 4_n . As a consequence, the Hilbert norm 
||j4|| 2 will be finite for n < |. We have thus succeeded in constructing an infinite class of 
zero modes, parametrized by the exponent 2 < n < I of eq. ([l9]) . whose corresponding 
gauge fields have finite Hilbert norm |[^4|| 2 . 



4 The case d = 4 

As last example, let us now face the case d = 4. As before, we adopt spherical coordinates 
(r,ct,9,ip). As starting configuration fulfilling the Landau gauge, we take 

Al = £ CF 4^W , d^Al = , /x = 1, 2, 3, 4 , (51) 

where g(r) depends only on the radial coordinate and e^vpa stands for the Levi-Civita 
antisymmetric tensor. Therefore, choosing the zero mode u a (r, a, 9, <p) in the usual way, 
namely 

u a (r,a,9,p) = (w (1) ,w (2) ,0) , (52) 
it follows that (u;W,u/ 2 )) obey the differential equations 

dV 1 ) + g(r)e 3flu4Xfl d u ojM = , 

8 2 ojW - g(r)e 3 ^x^u;W = Q . (53) 
Once more, one recognizes the appearance of the angular momentum in eqs. (|53p . i.e. 



d 

zzuvAX^dv = xd y - yd x = — . (54) 

Eqs. ([53|) reduce thus to 

a yi) + 9 ( r )i w p) = 0) 

dip 

5 y2)_ ff(r )i w (i) = . (55) 
dtp 

Generalizing the same procedure adopted in the previous cases d = 2,3, we parametrize 
the zero mode oj a as: 

u a (r, a, 9, if) = (<r(r) sin 9 sin a sin ip, a(r) sin 9 sin a cos <p, 0) . (56) 
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Equations ([55]) become a single equation, given by 



o"(r) + -a'(r) - \a(r) = g(r)a{r) . (57) 



Employing the change of variables c(r) = rip(r), eq. (|57|) yields 



g(r) = (<(r) + ^'M ) . (58) 



As usual, we set 



ri + a 

where q is a positive exponent and a a positive constant. From (|58|) . for the function g(r) 
we can write 

3^ = {rq l a) 2 (^~V ~ 4?) " ar«" V + 4g)) . (60) 
For the norm of the zero mode, we have 

,2 / ■ 



llwll oc / dr- -7T , (61) 

which is finite for q > 3. Expression (|59p provides thus an infinite class of normalizable 
zero modes. 

Concerning now the function g(r) in eq. (|60p . it behaves as g{r) ~ 1/r 2 for r — > oo, giving 
rise to a gauge field ^4 ~ 1/r. These configurations do not have finite norm ||^4|| 2 , although 
yield a finite action. In fact, in this case, for the field strength one has 

Fij = 2e ajii g{r) - e aik A x k x j ^— + e ajk4 x k x 1 ^— + e abc e bi k4 e cjl4 x k x l g(r) 2 , (62) 
which, upon contraction of the indices, becomes 

Ffj = e ajU (2g(r) + rg'{r)) + e ijki (g(rf + 9 -^)x a x k . (63) 
Therefore, for the Yang-Mills action we obtain 

J tfxFfjF* = 2vr 2 / drr 3 [6 (2g{r) + rg'(r)f - 4r 2 (2g(r) + rg'(r)) + g{rf 

+2r 4 (^+ 5 (r) 2 )] , (64) 
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which is easily seen to be finite, as g(r) ~ 1/r 2 . 



Let us end this section by showing that, as in the case d = 3, also in d = 4 it turns 
out to be possible to find an infinite class of zero modes exhibiting finite Hilbert norm 
P|| 2 . To that purpose, we first remark that, for the particular case q = 4, expression ([60]) 
becomes 

so that g ~ 1/r 6 for r — > oo. As such, the corresponding gauge field behaves as ~ 1/r 5 , 
exhibiting thus finite norm ||^4|| 2 . As done in the previous case d = 3, we consider now the 
generalized expression 

Mr) = —. , (66) 

where n is a positive exponent smaller than 4, n < 4, and /3, a are positive constants. From 
eq. ([58j) . the corresponding expression for g(r) is found to be 

g (r) = — — {(3 2 r 2n ~ 2 (2n 2 - 5n) + /3r n+2 (12n - n 2 - 32) 

( t I py I ) 

+r n - 2 fia(f3n - 5n - pn 2 ) - 32ar 2 } . (67) 

We see that, as in the case d = 3, we must constrain n > 2 in order to avoid the appearance 
of a singularity at the origin. Asymptotically, we have g(r) ~ l/r 6_n and thus Ai ~ l/r 5_n . 
As a consequence, the Hilbert norm ||vl|| 2 will be finite for n < 3. Thus, also in d = 4, we 
have been able to construct an infinite class of zero modes, parametrized by the exponent 
2 < n < 3 of eqs. (|66p . ([67p . whose corresponding gauge fields have finite Hilbert norm ||j4|| 2 . 



5 Conclusion 

In this work we have pursued the investigation of the issue of the Gribov copies. In partic- 
ular, we have been able to provide examples of infinite classes of normalizable zero modes 
of the Faddeev-Popov operator in the Landau gauge, in d = 2, 3, 4 dimensions, and for the 
gauge group SU(2). 

Two boundary conditions for the behavior of the gauge field at infinity have been consid- 
ered, namely: the requirement of finite Yang-Mills action, i.e. Sym = 4^7 / d A xF^F^ < 00, 
and that of finite Hilbert norm, i.e. \\A\\ 2 < 00. The requirement of finite Hilbert norm 
is stronger than that of finite Yang-Mills action. Gauge configurations leading to a finite 
norm automatically give rise to a finite action, the converse being not true. Nevertheless, 
in both cases, we have been able to provide in a rather simple way explicit examples of 
infinite classes of normalizable zero modes with the desired behavior for the gauge field at 
infinity. Here, Henyey's set up [3] has been proven to be particularly useful. 



12 



Let us conclude with a few remarks. 



• the examples of the zero modes and related gauge field configurations which we have 
been able to construct in a rather easy way give a simple confirmation of the relevance 
and, at the same time, of the complexity of the issue of the Gribov copies for a correct 
quantization of Yang-Mills theories. 

• it is worth to point out that the Gribov region f2 

Q = { A1 ; diAf = , M ab = - (d 2 5 ab + e acb A^ > } , (68) 

can be introduced without any specific reference to a particular behavior of the gauge 
field at infinity. It is apparent in fact that the region Q can be introduced in a way 
which is perfectly compatible with both boundary conditions which we have taken 
into account. In this sense, the restriction to that region, as implemented by the 
Gribov-Zwanziger action [HO US], has to be regarded as an important step towards 
a correct quantization of Yang-Mills theories. Certainly, the restriction to the region 
f£ enables one to eliminate all Gribov copies related to zero modes of the Faddeev- 
Popov operator. 

• As we have already underlined, several properties of the Gribov region fi, eq. ()68p . 
have been established [22], [231 123] , namely 

— i) the region Q is convex and bounded in all directions in field space, 

— %%) every gauge orbit crosses f2 at least once. 

In particular, property ii) provides a support to the original Gribov proposal of 
restricting the domain of integration in the functional integral to the region This 
restriction is achieved by the introduction of the Gribov-Zwanziger action [T| 114} [T5]. 
which can be summarized by the following formula 

/ [DA] 5{dA) det{M) e~ SYM = f [DA] 8{8A) det{M) e -( s ™+s„-/ <*W(iV 2 -i)) j 
Jn J 

(69) 

where Sh is called horizon function and is given by the nonlocal expression 

Sh = 7 4 J d A x h(x) , 

h(x) = g 2 J d 4 yf bal A^x)(M- 1 ) lm (x,y)f bkm A>;(y) . (70) 

The massive parameter 7 is the Gribov parameter, defined in a self-consistent way 
through the gap equation [14] [15] 

(h{x)) = 4 7 4 (iV 2 - 1) . (71) 
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The parameter 7 encodes the presence of the Gribov horizon, i.e. it originates from 
the restriction of the domain of integration in the functional integral to the region £1. 
As such, 7 can be seen as arising from the elimination in the functional measure of 
the field configurations which lie outside of the region f2, and which are necessarily 
Gribov copies of configurations belonging to £1. As shown in [14\ I15| . the nonlocal 
horizon function implementing the restriction to f2 can be cast in a local form through 
the introduction of a suitable set of localizing fields. Remarkably, the resulting local 
action turns out to be renormalizable to all orders. Recently, a refined version of the 
Gribov-Zwanziger action has been worked out |16 [ \17 \ \TE\. The predictions obtained 
from the Refined Gribov-Zwanziger action turns out to be in very good agreement 
with the lattice data on the gluon and ghost propagators, see [19j [20] and references 
therein, as well as on the spectrum of the lightest glueballs [21]. These preliminary 
results can be taken as evidence of the important role played by the Gribov issue 
towards a satisfactory understanding of the infrared behavior of confining Yang- Mills 
theories. 

• Finally, we should touch on the fact that the restriction to the Gribov region is 
not yet sufficient to eliminate all possible copies. Additional copies are still present 
inside the region fi. As already mentioned, there exists a smaller region, contained 
inside the Gribov region, which is free from Gribov copies. This region is known as 
the Fundamental Modular Region [221 123 1 [24"]. It is known that part of the boundary 
of the Gribov region is common to the Fundamental Modular Region. In princi- 
ple, the restriction of the domain of integration to the Fundamental Modular Region 
should provide a quantization formula for Yang-Mills theories free from any Gribov 
ambiguities. Unfortunately, unlike the case of the Gribov region Q, a local and renor- 
malizable action implementing the restriction to the Fundamental Modular Region 
is not available. In view of these remarks, it is safe to state that the possible physical 
relevance of the field configurations which we have investigated in this paper could 
be asserted only after a detailed study of their relationship with the Fundamental 
Modular Region, a rather difficult task which is out of the aim of the present work 
which attempts at providing examples of the zero modes of the Faddeev-Popov op- 
erator. For a better idea of the relevance of this issue we remind the reader to the 
lattice investigations given in [25j [26] . 
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Appendix 



For the benefit of the reader, we remind here some basic expressions in polar coordinates 
in d = 2, 3, 4 dimensions. 

d=2 

• polar coordinates 



x = r cos (p 

y = r sin ip (72) 



where ip ranges over [0, 2tt] 
• orthonormal basis of unit vectors 



e r = i cos ip + j sin ip 

e v = —i sin ip + j cos ip (73) 



differential operators 
gradient 



Laplacian 

Fi 2 f = V7 2 f = 

r dr \ dr J r 2 dip 2 



d=3 

• the spherical coordinates are given by 



x = r cos 93 sin i 
y = r sin 93 sin < 
z = r cos 6* 



where ip ranges over [0, 27r] and ranges over [0, ir] 



(76) 
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basis of orthonormal vectors 



e r = i sin 9 cos ip + j sin 9 sin ip + k cos 9 
eg = i cos 6* cos ip + j cos sin ip — k sin 
e<p = — « sin ip + j cos <p 



differential operators 
gradient 



* f .df ldf 1 df 
v / = e r— + efli^ + e v 



dr r 89 r sin 6 dp 
Laplacian 

~ 2f w2f 19 / 2 0A , i a ( . an i a 2 / 



r 2 <9r \ dr J r 2 sin 9 89 \ 89 J r 2 sin 2 9 dtp' 
[ 

in four dimensions, for the spherical coordinates we have 

x = r cos p sin 9 sin a 
y = r sin <p sin sin a 
z = r cos # sin a 
t = r cos a 

where ip ranges over [0, 2tt], while (9, a) range over [0, ir] 
basis of orthonormal vectors 



e r = i sin cos <p sin a + j sin sin ip sin a + k cos sin a + t cos a 

eg = i cos cos <p + j cos sin p — k sin 

e^, = — i sin 99 + j cos 

e a = i sin cos <p cos a + j sin 6* sin ip cos a + k cos cos a — i sin a 



differential operators 
gradient 



V/ = e r — + e Q -— + eg — : — + e v — : . 

or r aa r sin aw r sin a sin a<^ 
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Laplacian 



d 2 f = v 2 / = I- fr 3<9/> ) I 1 ^ fsin 2 | 1 f 1 9 (sine 9 !) + 1 ^ 

r 3 dr \ dr J r 2 sin 2 a da \ da J r 2 sin 2 a \sm8 89 \ 86 J sin 2 9 dip 2 

(83) 
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